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Abstract. After a brief discussion of effective field theory applied to nuclear clusters, I concentrate on the 
inclusion of two particular aspects, namely, narrow resonances and electromagnetic interactions. As examples of 
applications, I present the details of our studies on alpha-alpha and proton-alpha scattering. 



1 Introduction 

Two-particle scattering at low enough energies are insen- 
sitive to the details of the interparticle interaction, pro- 
vided that the latter has a limited range of action. The in- 
tuitive argument behind this general property is that mo- 
menta whose wavelengths are much larger than the inter- 
action range cannot resolve its finer details. The elastic am- 
plitude incorporates unitarity constraints and an universal 
low-energy expansion, a.k.a. the effective range expansion 
(ERE), whose coefficients encode the low-energy behavior 
of the forces for a particular system. This simple observa- 
tion, made by Bethe and others 1 1 1 more than fifty years 
ago, became a benchmark for theorists trying to model the 
dynamics of the interaction. 

The same line of reasoning is shared by effective field 
theories (EFT), where only the relevant low-energy de- 
grees of freedom are explicitly taken into account. Modes 
that are active only at high energies are "frozen" into the 
low-energy constants. In the particular case of EFT with 
short-range interactions, the two-body scattering amplitude 
is equivalent to a low-energy expansion of the ERE am- 
plitude (2). However, the EFT approach has extra features 
over the effective range theory. First, the systematic expan- 
sion of the amplitude allows for a reliable and unbiased 
estimate of the theoretical uncertainties. Second, symme- 
tries of the problem or the way they are broken can be im- 
plemented in a straightforward way (for instance, electro- 
magnetic and weak currents, isospin, parity). And third, the 
ability to extend the formaUsm to systems with more than 
two particles. 

One example of the second feature illustrates a differ- 
ence between the EFT and ERE approaches. In the pro- 
cess np — » dy, the Mly transition contains at NLO a four- 
nucleon-photon operator that cannot be generated by the 
ERE approach |3|. The associated coupling can be fixed 
from cold neutron capture rate at incident neutron speed of 
2200 m/s. Rupak improved the calculation of Mly and Ely 
transitions to N^LO and N^LO, respectively, lowering the 
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theoretical uncertainty in the np ^ dy reaction to about 
1% for center-of-mass energies up to 1 MeV L4J- 

When extended to three or more particles with large 
scattering length(s), EFT provides new insights and ex- 
plains certain universal properties that are independent of 
the interaction details. For instance, the related Thomas 
and Efimov effects |5 1 are intimately connected to the be- 
havior of the leading three-body counterterm under varia- 
tions of the momentum cutoff in the dynamical equations. 
Such counterterm is necessary to properly renormalize the 
theory, as well as to guarantee a unique solution in the 
limit when the cutoff is taken to infinity f6'71. It exhibits 
a limit cycle behavior, that is, a log-periodic dependence 
with the cutoff, that explains the geometrically-separated 
bound states in the Efimov spectrum. For nuclear systems 
with three and four nucleons, the theory provides a con- 
vincing explanation to certain correlations, like the almost 
linear dependence of the spin-doublet neutron-deuteron scat- 
tering length with the triton binding energy II8I9I or simi- 
lar dependence of the latter with the alpha particle binding 
energy ifTOl , known respectively as the Phillips and Tjon 
lines. 

Due to the universal character, effective field theories 
with short-range interactions have been successfully ap- 
plied to distinct areas of physics — atomic, particle, and 
few hadron systems. For a more complete overview, see 
Refs. M6I1 II . Here I will concentrate on applications to 
nuclei that behave as systems of loosely bound clusters, 
which also include the nowadays popular halo nuclei. In 
Section |2]the general features of the so-called halo/cluster 
EFT are presented in certain detail, with emphasis on as- 
pects quite often faced when dealing with nuclear clusters. 
Sections |3] and m present applications of halo/cluster EFT 
to alpha-alpha (aa) and proton-alpha {pa) scattering, and 
in Section|5]are the concluding remarks. 



2 EFT for nuclear clusters 

The relevant degrees of freedom in halo/cluster EFT are 
the structureless, weakly bound objects represented by re- 
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spective field operators. In the case of halo nuclei, these 
are the core nucleus and the valence nucleons. Effects like 
nucleon excitations inside the core, pion or nucleon ex- 
changes, take place at energies well above the binding en- 
ergy that holds the clusters together. The former has an 
associated momentum scale M/„ of the order of the inverse 
of the interaction range, while in the latter case the mo- 
mentum scale Mio <sc M;,, is inversely proportional to the 
size of the halo system. As explained in the following, this 
separation of scales is, from the EFT point of view, the 
outcome of a fine-tune in the coupling constants, similar to 
the one responsible for generating a shallow bound state 
in the spin-triplet nucleon-nucleon channel. In fact, due 
to its large extension compared to the typical range of the 
proton-neutron interaction, the deuteron can be seen as the 
simplest halo nucleus, with a loosely bound neutron sur- 
rounding a proton core. 

In order to understand the formalism, let us consider a 
system of identical bosons, with mass m^, represented by 
a field (p and interacting via a pairwise S -wave short-range 
force. The most general Lagrangian respecting the relevant 
symmetries of the system (parity, total angular momentum, 
approximate isospin and non-relativistic Galilean invari- 
ance) is given by 
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where we introduce an auxiliary (dimeron) field d, with 
"residual mass" A, carrying the quantum numbers of two 
bosons in 5 -wave and coupling with their fields through 
the coupling constant g. The dots stand for higher order 
terms in a derivative expansion. The building blocks are 
the boson and dimeron propagators, the LO coupling g 
between the dimeron and two boson fields, and higher- 
order couplings that can be rewritten, via field redefini- 
tions, as powers of the dimeron kinetic energy. The boson 
and dimeron propagators read 



iS a(qo; q) 
iD(qo; q) 
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where the precise form of the latter depends on the magni- 
tude of each term in the denominator The power counting 
specifies the amount of fine-tuning in the low-energy cou- 
plings, responsible for assigning a well-defined order to ev- 
ery possible Feynman diagram that contributes to the scat- 
tering amplitude. In a natural scenario, where fine-tuning 
is absent, the effective couplings in Eq.([TJ depend only on 
the high-energy scale M/,, of the theory. That corresponds 
to zf ~ Mf^./nia and ~ Mhi/m^. The bare dimeron prop- 
agator (Is) is static at leading order (LO), 
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The LO tree diagram is proportional to aq = niag^/AnA ~ 
0(1 1 M hi), where qq is the boson-boson scattering length. 



The one-loop graph receives two additional factors from 
the coupling constant g, one dimeron propagator, one-loop 
integration measure dq^d^q ~ k^/nia, and two interme- 
diate boson propagators 5^ ~ {nia/k^)^. The net result is 
a suppression of k/Mi,j compared to the tree graph, with 
k the typical low-energy momentum under consideration. 
The next order, k^/Mf^. smaller than the tree level, com- 
prises two-loop contributions proportional to (aok)^ and 
one kinetic term insertion to the dimeron propagator pro- 
portional to a^rok^, where ro = A-n/mlg^ ~ l/M/„- is the 
boson-boson effective range. The amplitude then amounts 
to a simple Taylor expansion in powers of the expansion 
parameter k/M^i- 

However, a natural perturbative treatment is unable to 
account for a plethora of nuclear phenomena, where bound 
states and resonances are the rule rather than exceptions. 
That demands a certain amount of fine-tuning in the pa- 
rameters of the effective Lagrangian. The physics of shal- 
low S -wave bound states can be described by fine-tuning 
the parameter A in Eq.([T]i, responsible for driving the size 
of the scattering length a^. Setting A ~ Mi,iMio/ma, leads 
to a large oo ~ 1 /Mi,,, with Mio the new momentum scale 
associated to the fine-tuning. This situation is analogous 
to the nucleon-nucleon (A^A^) case |6 7|. A consequence of 
this new scale is that higher-loop graphs built up from the 
LO boson-boson-dimeron coupling g are no longer sup- 
pressed when k ~ M/,,, but contribute at the same order as 
the lower-loop and tree graphs. The resummation of this 
class of diagrams leads to a LO amplitude in the form of 
a truncated effective range expansion up to its first coef- 
ficient. The scattering length then becomes the only rele- 
vant scale and (for positive values) predicts the existence 
of a shallow bound state Eb ~ lliniaa^^. The renormal- 
ization group (RG) analysis of the Schrodinger equation 
demonstrates that two-body systems with large scattering 
length consist of strongly interacting systems with an RG- 
flow towards a non-trivial fixed point |12|, in contrast to 
the natural case of weakly interacting particles, where the 
RG flows towards a stable, trivial fixed point. 

When extended to three-body systems, this EFT pro- 
vides new insights on some universal features, as men- 
tioned in the introduction. An investigation searching for 
these universal effects in two-neutron halos was performed 
in Ref. [ 13 1, where it was found the ^"C nuclei as possible 
candidate to have an Efimov-like spectrum (see also [il4|). 
Non-universal (effective range) corrections have recently 
being incorporated to this study ifTSl . 

There is an additional phenomenum quite often present 
in nuclear clusters, namely, the existence of low-energy 
resonances in the continuum. One could then expect the 
need for extra amount of fine-tuning. This is the case for 
neutron-alpha {na) scattering II16I17I . where the presence 
of a narrow resonance in the f 3/2 channel around E„ « 1 
MeV requires adjustment of the coupling constants of the 
theory to enhance the P3/2 amplitude a few orders of mag- 
nitude away from the natural assumption. Another exam- 
ple is the aa system [18], which contains an 5 -wave reso- 
nance about 0. 1 MeV above the elastic threshold. The aa 
scattering amplitude can be described by the Lagrangian 
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([TJ, but assuming A ~ Mf^Jiria- This is necessary to gen- 
erate an even larger aa scattering length, ao ~ Mhi/Mf^. 
These scalings will be discussed in detail in the next sec- 
tion. With this power-counting the dimeron propagator is 
no longer static, as the kinetic and residual mass terms are 
now of comparable order In momentum space it has the 
form of Eq.(|3]l. The presence of the kinetic term in the 
dimeron propagator in practice resums effective range ro 
corrections to all orders. This resummation is necessary to 
reproduce a narrow resonance at low energv lll8ll6ll7l . 



+ 



+ 



Fig. 1. Graphic representation of Tcs, as sum of graphs with 
multiple insertions of the bare dimeron propagator (double line) 
and the "Coulomb bubble loop". The latter contains intermediate 
Coulomb photons resummed to all orders (shaded ellipse). 



For the case of charged bosons, electromagnetic inter- 
actions can be incorporated into Eq.([T]) via the usual min- 
imal substitution. Among them. Coulomb photons are the 
dominant ones at low energies 1 18|. Coulomb interactions 
were formulated in the EFT framework by Kong and Ravn- 
dal ||T9l for the two-proton system, and can be extended in 
a straightforward way to include resonances 1 18 1. The idea 
relies on the two-potential formalism, where the amplitude 
is separated in a term that contains only the two-particle 
pure Coulomb scattering Tc, plus the Coulomb-modified 
strong amplitude Tcs ■ The latter is diagramatically illus- 
trated by Fig.[T](see caption) and has the form of a geomet- 
ric series. Like in the proton-proton case lfT9l . the power 
counting for narrow resonances demands its resummation. 
Up to next-to-leading order (NLO) one gets for our two- 
boson example, with reduced mass fi and charge Zq,, 
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where 
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is the inverse of the Bohr radius, and Pq, the shape pa- 
rameter. The strength of Coulomb photons is driven by the 
Sommerfeld parameter 
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and the Coulomb-modified strong amplitude is modulated 
by the Sommerfeld factor. 



ci 



(7) 



which quantifies the probability of finding the two bosons 
at the same point in space. The phase of pure Coulomb 
relative to free scattering, for a given partial wave L, reads 



(Tl = arg r{L + 1 + irf) 
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The last Coulomb ingredient in Eq.(|4]i is the function 
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and explicitly shows the complicated analytic structure of 
Coulomb scattering at low energies. Nevertheless, apart 
from the Sommerfeld factor and the S -wave Coulomb phase, 
one notices that the effect of dressing the strong amplitude 
with non-perturbative Coulomb photons amounts to, apart 
from a multiplicative factor in the numerator, replacing the 
unitarity term -ik by -2kc H(rf) in the denominator of the 
amplitude. Similar expression can be derived for P-wave 
interactions and will be presented in Section |4] 



3 aa scattering 

3.1 power counting: pinning the scales down 

In order to better understand the scalings of the EFT pa- 
rameters in the presence of a resonance, it is helpful to 
provide some numbers for the relevant physical quantities. 
It is well-known that the aa system is dominated by 5- 
wave at low energies |20|, having a very narrow resonance 
at Er ^ 92 keV and width Fr x 6 MeV. To this reso- 
nance is commonly assigned the ^Be ground state. The aa 
scattering length aq ~ 2000 fm is about three orders of 
magnitude larger than the a matter radius, suggesting the 
large amount of fine-tuning discussed in the last section. In 
contrast, the effective range ro ~ 1 fm and shape parameter 
Pq ~ -1 .Ifnr' obey the expected natural scalings of 1 /M/„- 
and ^/M^., respectively. It seems therefore reasonable to 
start with the scalings A ~ M^^/jj and g^/2n ~ Mi,i/fi^. 
With the evaluation of Feynman diagrams and after the re- 
quired resummation described in the previous section one 
finds 



Tcs---Cy-" 



27r^ 



Mr 



^k^ - 2kcH{iT) 



(10) 

where /I**' is the only parameter in the Lagrangian, up to 
NLO, that is renormalized by the Coulomb loops. 
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with D the number of space-time dimensions, k the renor- 
malization scale parameter, and Ce - 0.577 . . . the Euler- 
Mascheroni constant. The amplitude matches with the ex- 
pansion of the ERE amplitude (|4]l, from where one iden- 
tifies the ERE parameters in terms of the EFT couplings. 
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One then finds gq ~ Mi,j/Mf^^. The low-energy scale Mio ~ 
yjmaER a: 20 MeV is nearly seven times smaller than a 
high momentum scale associated to either the pion mass 
or the excitation energy of the alpha particle, M;„- ~ m„ ~ 
-^JniaE* X 140 MeV. Based on these numbers, one could 
expect convergent results for observables at laboratory en- 
ergies up to 3 MeV. 

The Bohr momentum kc provides the scale of Coulomb 
interactions. Due to the large aa reduced mass in its def- 
inition (|5]l, kc turns out to be numerically of the order of 
Mhi- One can therefore expects large and important elec- 
tromagnetic contributions. Nevertheless, it is interesting to 
discuss the other limit kc — > 0, when Coulomb is turned 
off. In this case one has 2kcH{r]) — > ik and the denomina- 
tor of the LO amplitude reads 



-^LO °^ ~l/flo + f()k^ 12 — ik . 



(12) 



For momenta k ~ Mio the first two terms are suppressed 
by MiolMhi compared to the last one. Therefore, all that 
is left at LO is the unitary term l/(-ik). In this limit, the 
^Be system exhibits non-relativistic conformal invariance 
pij and the corresponding three-body system, '^C, ac- 
quires an exact Efimov spectrum 1.6.8] . If it was possi- 
ble to somehow shield the charges of the a particles, the 
above limit strongly suggests that two and three charg- 
less a particles would probably be the best nuclear systems 
to observe few-body universality. The above scenario cer- 
tainly changes for a physical value of kc, due to the break- 
ing of conformal invariance by the 1/r Coulomb force. 
Nevertheless, the fact that the ground state of ^'Be and the 
Hoyle state in ^^C remain very close to the threshold into 
a-particles suggests that this conformal picture is not far 
from the real case. 

The breaking of scale invariance by the Coulomb in- 
teraction introduces the scale kc in the propagation of two 
charged particles. As we have seen in Eq. (|4]i, the unitarity 
term is modified. The balance between strong-interaction 
terms and Coulomb-modified propagation now depends on 
both the strong-interaction scale Mio and kc- While the for- 
mer is kR = s/niaER ~ 20 MeV, the latter is kc = aemZ^fi ~ 
60 MeV, and therefore a relative strength of kc/kn ~ 3. 
For momenta k ~ k^, we are clearly in the deep non- 
perturbative Coulomb region. The Sommerfeld parameter 
T] reaches large values and the function IkcHirf) is very 
diff'erent from the usual unitarity term ik. Instead of ham- 
pering any simplification, it actually allows a low-energy 
expansion of the Coulomb function H. Using Stirling's se- 
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and tl/{z) = (d/dz)\nr(z) in Eq. Q gives 
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The unitarity term is thus replaced by 2kcH(ri) ~ k^/6kc at 
LO. This is now a factor k/6kc smaller in magnitude than 



the unitarity term in the absence of Coulomb, and compa- 
rable to the eff'ective-range term coming from the dimeron 
kinetic term. One can grasp it automatically if one consid- 
ers 3kc ~ Mhi, as it appears to be the case numerically. 



3.2 resonance pole expansion 

Before incorporating the simplifications on the function 
Hirff, one should take a careful look at the EFT amplitude 
up to NLO. First, an isolated remark. Since we assume that 
!P() is a parametrically small correction, one is allowed to 
reshuffle the perturbative series and resum Vq to all orders, 
thus recovering the ERE formula. 
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Second, one notices that equation (HJi is valid for generic 
momenta k ~ Mi,,. However, it fails in the immediate prox- 
imity of k^. This situation is familiar from the neutron- 
alpha case IfTTl . The power counting works for k ~ Mio 
except in the narrow region \k - kfi\ = 0(Mf^/Mhi) where 
the LO denominator approaches zero and a resummation 
of the NLO term, here associated with the shape parameter, 
is required. As one gets closer to the resonance momentum 
k]i, higher-order terms in the ERE are kinematically fine- 
tuned as well. This happens because the imaginary part of 
the denominator is exponentially suppressed by a factor 
exp(-2;T77R) ~ 10"^ and the real part is allowed to be ar- 
bitrarily small. Nevertheless, this kinematical fine-tuning 
is not a conceptual problem. From the EFT point of view, 
each new fine-tuning can be accommodated by reshuffling 
the series and redefining the pole position. Such a proce- 
dure works fine with a small number of kinematical fine- 
tunings, but is not practical in the aa system. A better al- 
ternative is to perform an expansion around the resonance 
pole position, starting from the resummed (ERE) ampli- 
tude. The situation here is nearly identical to the A^A^ sys- 
tem, where one can choose to expand the amplitude around 
the bound-state pole Il22l rather than around zero energy. 

A great simplification is achieved from the fact that the 
resonance lies in the deep Coulomb regime, where Eq. (fT4l) 
provides an accurate representation of H up to the preci- 
sion we are considering. The real terms shown in Eq. ( fT4b 
then become an expansion in powers of ~ {k/3kc)^ = 
0(k^/Mj^.). Given the asymptotic feature of the Stirling's 
series, at some point the remaining terms can no longer be 
expanded; at that point the remainder should be treated ex- 
actly. In lowest orders, however, we can use the successive 
terms shown in Eq. (fT4l i: numerically, the terms up to ;/ '* 
work to better than 3% for Eu^b - 3 MeV. 

The expansion ( fT4l i not only simplifies a complicated 
function of kc/k, but also makes the physics around the 
resonance more transparent. Since the "size" of the res- 
onance, l/ku, is much larger than the Bohr radius 1/kc, 
the Coulomb interaction is eff'ectively short ranged, and 
the real part of H resembles the ERE expansion. In the 
amplitude Tcs the different strong and Coulomb coeffi- 
cients proportional to a common power of k can be grouped 



19 International lUPAP Conference on Few-Body Problems in Physics 



together, where one observes that they have comparable 
sizes. We therefore define 



1 



1 



and so on. Up to NLO we rewrite Tcs as 
In 
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From this expression one sees directly that k/t ~ Mu,, with 
corrections of 0{Mf^/Mf^.). Note that we keep the exact 
form of the imaginary term in Eq. (fTTl i: even though it is 
negligible at A: ~ k^, it has an important exponential de- 
pendence on the energy responsible for keeping the phase 
shifts real in the elastic regime. 

When fl() < and ro < l/3kc, we have kj, > and 
the two poles of Eq. ( fTTI l are located in the lower half of 
the complex-momentum plane very near the real axis, as 
expected for a narrow resonance. The amplitude Tcs can 
be written in terms of the resonance energy Er = kg/2fi 
and the resonance width r(E) a£] 

Tcs = ^ ^ ■ (19) 
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One finds that 
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See also Ref.a23i. 



The width is very small because of the large value of 2nkclkR 
in the exponential. 

In the form of Eqs. ( fT9] l and ( l20t we can keep E^ and 
r{ER) fixed at each order in the expansion. Note that these 
equations do not change to this order if one makes a dif- 
ferent choice — e.g., (k^ - k^)^ instead of k'^ - k'^ — for the 
form of the !Po term in Eq. (fTTl l. The behavior of the phase 
shift around Er is guaranteed to be of a resonant type, as it 
automatically satisfies the constraints 
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(resonance energy) , (22) 



(resonance width) . (23) 



3.3 confronting the data 

Unfortunately aa scattering data at low energies are not 
abundant. Nevertheless, all the existing measurements at 
Elab up to 5 MeV show that it is dominated by the S wave, 
thanks to the presence of the (J'', I) = (0^, 0) resonance im- 
mediately above threshold. Early determinations of the 0^ 
resonance energy were performed in reactions like " B H- 
p ^ 2a + a (see [20] and references therein). Later mea- 
surements of the scattering of "^He atoms on "^He^ ions, es- 
pecially projected to scan the resonance energy region and 
supplemented with a detailed analysis |24 25| improved 
the determination of the resonance energy and width to 
their currently accepted values, E^^^ - 184.15 ± 0.07 keV 
and rf^g = 1 1.14 + 0.50 eV 125). The resonance CM mo- 



mentum is thus kji - ^fiEf^^g » 18.5 MeV. 

The obvious way to extract the ERE parameters from 
the 5 -wave phase shift is to fit (the cotangent of) the lat- 
ter to the ERE formula. However, fits to aa scattering data 
alone are not able to constrain well those parameters and, 
apart from the large uncertainties, don't seem to predict the 
existence of the 0^ resonance. Imposing that the phase shift 
crosses the value 7r/2 at the resonance energy improves the 
situation, though not dramatically. Following a previous 
suggestion 1261 . Ref. f2l \ used not only the available res- 
onance energy, but also the width |24| to constrain these 
parameters. The obtained values for aq - -1.76 x 10^ fm, 
ro = 1 .096 fm, and = - 1 .654 fm^* made use of the reso- 
nance energy and width from Ref |24|. Later we compare 
these numbers with the ones from our EFT fits. 

At energies below Ei^j^g — 3 MeV data were obtained, 
and a phase-shift analysis performed, by Ref. Il28ll . Values 
of the latter can be found in Table II of the review |20|. 
Since it is well-determined experimentally, and due to its 
relevance to the triple-alpha process, we use the 0^ reso- 
nance parameters from Ref. |25| as important constraints. 
This is in line with the EFT approach, where lower-energy 
observables have preference over higher-energy ones. The 
relationship among the EFT parameters and the resonance 
energy and width allows one to reduce the number of vari- 
ables to be adjusted at each order in the power counting. 
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Below, we also show the ERE from Ref. ||27| for orien- 
tation, and comment on the extremely large value of the 
scattering length aq, which suggests a large amount of fine- 
tuning in the parameters of the underlying theory away 
from the naturalness assumption. 

In the power counting that we discussed for the aa sys- 
tem, the amplitude Tcs for generic momenta is given up to 
NLO by Eq. (fTTI ). As demonstrated in the previous subsec- 
tion, this expression combines the deep-non-perturbative 
Coulomb approximation (fT4l l for the function H with the 
expansion around the resonance pole, thus preventing the 
need for multiple kinematical fine-tunings. In LO, the two 
parameters and ro can be obtained from the constraints 
(fTSl l and dirt . At NLO, a fit to scattering data is needed to 
determine Pq. 




12 3 

Elab [MeV] 

Fig. 2. EFT results for aa scattering at LO (dotted) and NLO 
(solid), compared against the data. Top panel: phase shift S^. Bot- 
tom panel: K(r]) = C^(cot^^ - i)/2j] + H{r)). 

Figure |2] shows our results with the resonance position 
and width constraints, compared to the available 5 -wave 
phase shifts below Eiab - 3 MeV. In the region above the 
resonance, where scattering data are shown, the LO curve 
is a prediction, which is consistent with the first few points 
but then moves away from the data. The NLO curve has P() 
as an extra parameter, which is determined from a global 
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X^-^^ to scattering data shown. As expected from a conver- 
gent expansion, the description of the low-energy data im- 
proves with increasing order. At about 3 MeV and above, 
higher-order contributions are expected to be important, 
as suggested by the discussion on the relevant scales and 
manifest in the growing difference between the NLO curve 
and both LO curve and data points. Also shown are results 
from a fit using the conventional ERE formula, Eq. ( fTsl l. 
in order to stress the differences between this and our EFT 
approach. 

Table [1] shows the values of ERE parameters used to 
produce the curves of Fig. |2l and compare them with the 
values 1 27 1 obtained from effective range theory. At LO, 
our values for oq and ro are consistent with the ones from 
Ref. f27L remarkably ro. The NLO corrections, however, 
spoil this initial LO agreement. The reason for this devia- 
tion could be due to the way the width constraint was ob- 
tained in Ref. f2T|. Its Eq. (4) reads 

dh 1 n 1 / „ T\ 

dk^ l^r{E]i) e2Mr/*/< - 1 4kc ^ ' 

where h{r\) = Re[H{ri)]. Following this reference's pro- 
cedure, we reproduced the quoted value of the width (6.4 
eV) only when dh(rij^)/dk^ was approximated by l/12k^ 
[see Eq. (fT4li1. However, that is equivalent to neglect the 
electromagnetic term of Pq in Eq. (|2TI) . which is incon- 
sistent since the strong piece contributes at the same order. 
Neglecting this term explains the agreement at LO and dis- 
agreement at NLO between our and Ref. lIZTl numbers for 
ro. With a smaller ro, and therefore larger (negative) ro, 
one can also understand why Ref. [27 1 obtains a smaller 
flo, as the product rooo is inversely proportional to the res- 
onance momentum squared. By repeating the Ref. |27|'s 
procedure including the width constraint consistently, we 
obtained essentially the same values as in our EFT fits. 
This updated ERE fit is also shown in Table [T] 





ao (10^ fm) 


ro (fm) 


n (fmO 


LO 


-1.80 


1.083 




NLO 


- 1.92 ± 0.09 


1.098 ±0.005 


-1.46 ±0.08 


ERE (our fit) 


- 1.92 ± 0.09 


1.099 ±0.005 


-1.62 ±0.08 


ERE |27| 


-1.65 ± 0.17 


1.084 ±0.011 


-1.76 ±0.22 



Table 1. ERE parameters extracted from EFT fits in the first two 
orders, compared with values from two ERE fits, our own and 
Ref. I27l's. 



Our fits reveal both effective range ro ~ 1/(180 MeV) 
and shape parameter ~ 1/(170 MeV)^ scaling with 
powers of M/„-, in agreement with our a priori estimate. 
The relative errors in oo and ro at LO are estimated to be 
of the order of the EFT parameter expansion, M/o/M/„ ~ 
1 /7 15%. At NLO, the main source of uncertainty comes 
from the precision of the most recent measurement of the 
resonance width [25 1, lying between 4-5%. The uncer- 
tainty in To, given by the;('^-fit, is of the same order. Note 
that the small relative error in ro compared to the one in 
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ro is due to the former being an order of magnitude larger 
than the latter, as we discuss in the next subsection. The 
NLO errors found here are a factor of two smaller than the 
ones obtained by Ref. |27|. 

One should stress that the accurate value of the reso- 
nance width, Fr - 5.57+0.25 eV, imposes tight constraints 
on our fits, through qq and ro. A significant improvement in 
our NLO fit and overall agreement with data is observed. 
But looking in detail, one sees that the theoretical curve is 
not able to cross the error bands of many scattering points 
below 3 Me V. This can be inferred from the ;^'^/datum ^ 4. 
In principle, a better agreement should be achieved by an 
N^LO calculation. However, that would introduce an ex- 
tra parameter that is mostly determined by the scattering 
data, and an expected agreement could mask any possi- 
ble inconsistencies between the phase shifts and the reso- 
nance parameters. The high NLO ;if^/datum suggests that 
they are not compatible with each other or, at least, one of 
them has overestimated precision. As a test we fitted both 
our NLO EFT and ERE expressions to the scattering data 
without the constraints from the resonance width. In these 
two cases, description of S -wave phase shifts is much bet- 
ter at the expense of an underpredicted resonance width, 
Fr = 4.9 + 0.6 eV with ERE and Fr = 2.87 + 0.23 eV with 
EFT. The ERE result is still consistent with the measured 
Fr thanks to its large error bar. In EFT, where lower-energy 
data have higher priority, the discrepancy is amplified. The 
problem is even more pronounced if the fit is performed 
using data up to 2.5 MeV instead of 3 MeV: the results 
Fr = 4.2 + 0.6 eV with ERE and Fr = 2.93 + 0.34 eV with 
EFT fall beyond the quoted experimental error bars. Oddly, 
this tendency continues as one lowers the upper limit in the 
fit. Reanalyses of the existing low-energy data or even new 
measurements seem necessary to resolve this discrepancy. 



3.4 fine-tuning puzzle 

A surprising feature of the aa system is the very large mag- 
nitude of flo, even if compared to the large scattering length 
observed in the nucleon-nucleon system. The latter case is 
thought to be the outcome of a fine-tuning in the QCD pa- 
rameters, giving rise to an anomalously low momentum 
scale. It is plausible to expect that this fine-tuning propa- 
gates to heavier systems. However, the enormous value of 
ao in the aa is suggestive of a more delicate tuning, with 
electromagnetic interactions playing a crucial role. 

To better understand the puzzle it is worth looking in 
details at Eq. ( fTTT i. The Coulomb loop contributions, pro- 
portional to the curly brackets, scale as M^^Jfi, given that 
g^lln ~ Mkil^p- and 2kc ~ M^. For the scaling ~ 
Mj^/fj. to hold the scale-dependent parameter /I (a-) must be 
of the same size of the Coulomb loops, and strongly can- 
cel each other to produce a result (M/o/M/„)^ ~ 100 times 
smaller This amount of fine-tuning is necessary to obtain 
a resonance at the right position. Similar cancellation is 
observed for the P-wave resonance in na |17|. However, 
in the present case the fine-tuning is caused by a delicate 
balance between the strong and electromagnetic forces. 



This fine-tuning scenario becomes even more enigmatic 
by considering the width. The tiny value of ^6 eV 
is intimately related to the resonance momentum, as can 
be immediately seen in Eq. ( ISTb . That is simply an effect 
of the Coulomb repulsion — the resonance trapped inside 
the envelope formed by the strong plus Coulomb potential 
has to tunnel through a larger potential barrier as its en- 
ergy gets smaller. The probability of tunneling is propor- 
tional to the resonance width, from where the smallness 
of Fr follows. Despite the tiny value, the width has an as- 
sociated scale of AnkcliJ-h) which is quite large. A natural 
estimative for this quantity would be ~ ^lill^^ implying 
that ro ~ i/Mhj. However, from the width constraint one 
finds ro ~ Mio/Mf^., leading to an ro that roughly cancels 
against l/3A:c with a remainder about 10% smaller From 
Eq. ( fTsT i one sees that the extra fine-tuning in ro is respon- 
sible for an extra increase in the scattering length by a fac- 
tor of Mi,i/Mio, becoming eff'ectively |aol ~ ^li/^h- 
remarkable, that if the strong forces generated an ro 11 % 
larger the ^Be ground state would be bound, with drastic 
consequences in the formation of elements in the universe 
(see also Ref.|29|). 



4 pa scattering 

The pa scattering at low energies is mostly dominated by 
the P3/2 partial wave, due to the presence of a resonance at 
proton energies Ep ^ 2 MeV. S -wave also gives an impor- 
tant contribution for Bp < 5 MeV, which is the low-energy 
region we are interested in (see below). In the P1/2 wave 
the phase shift varies smoothly up to Ep - 18 MeV, sug- 
gesting the existence of a very broad resonance for proton 
energies between 8 and 15 MeV |30|. However, to the or- 
der in the power counting and energies that we consider, it 
behaves as a typical higher-order effect. Higher wave phase 
shifts are less than a few fractions of degree below Ep ^ 6 
MeV, where D-waves start being relevant. Therefore, only 
two partial waves contribute to the low-energy EFT for the 
pa system, up to and including NLO. 

The Lagrangian for nucleon and alpha particles inter- 
acting via S 1 /2 and P3/2 partial waves is written as 1.31.16.171 



Na 



D2 




D2 




iDo + - — 
Irria 


<p + N'' 


iDo + 


N 



1 



-A 



0+ 



iDo + 



D2 



f'S^ ■ D(N<f>) + H.c. 

D2 



2(OTq, + OTfti) 

ND(f> - iDN)<p 



H.c. 



s'N4> + (p'N's 



iDo + 



iDo + 



2(ma + mff) 

D2 



2{ma + nij^) 



t. 



(25) 
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Table 2. Pj/t pa resonance parameters from the extended R- 
matrix analysis of Ref. 1321 . 



Table 3. S 1/2 and pa parameters extracted from the ERE fit 
of Ref. (33). 



k,. (MeV) 


ki (MeV) 


Er (MeV) 


r«/2 (MeV) 


51.1 


9.0 


1.69 


0.61 



where (p and are the alpha and nucleon fields with masses 
nia and m^, respectively, and 



(26) 



O0+ (fm) 


ro+ (fm) 


4.97 ±0.12 


1.295 ±0.082 



ai+ (fm^) 


n+ (fm-') 


^1+ (fm) 


-44.83 ±0.51 


-0.365 ±0.013 


-2.39 ±0.15 



in the P3/2 channel, and 



The auxiliary dimeron fields f and s couple to Na in P3/2 
and 5 1/2 waves, with leading-order coupling constants gi+ 
and go+, respectively. The S'/'s are 2x4 spin-transition ma- 
trices between J - 1/2 and J - 3/2 total angular momen- 
tum states. The sign variables ^o+, fi+ - ±1 in front of 
the kinetic term of the dimeron fields are adjusted to re- 
produce the signs of the respective effective ranges. The 
gauge-covariant derivative is defined in the usual way in 
terms of the photon field A^,, 



— dfi + ieZ 



1 + T3 



(27) 



where Z is the charge of the "particle" whose field the co- 
variant derivative acts on, T3 is the z-direction Pauli matrix 
acting in isospin space, and is the photon field. 

The power counting in essence consists of establishing 
how the size of the (renormalized) EFT coupling constants 
depend on the momentum scales present in the system one 
wishes to describe. The former have a direct relation to ob- 
servables, namely, the ERE parameters. In the pa system 
one is interested in describing the scattering region around 
the f 3/2 resonance located slk - kp ~ Mio, where 



kp — kj- iki 



(28) 



is the (complex) resonance momentum in terms of real 
quantities kr and A;,. Table |2] lists the resonance parame- 
ters extracted from the so-called extended /J-matrix analy- 
sis used in Ref. 1321 , which indicates that k^ ~ Mio ~ 50 



MeV and ki ~ MfjMhi ~ 10 MeV. The suppression of 
ki relative to k^ matches with narrow character of the f 3/2 
resonance. As in the aa case |18|, M/„ is set as the mo- 
mentum required to excite the a core or the pion mass, 
both around 140 MeV. Our expansion parameter is of the 
order of 1/3, allowing us to expect convergent results for 
proton energies up to 5 MeV. 

Table [3] shows the numerical values of the 51/2 and 
f 3/2 ERE parameters taken from Ref. Il33l . As we can no- 
tice, while !Pi+/4 ~ ro+/2 ~ l/M/,,- are in agreement with 
the natural assumption, the quantities 00+ ~ ai+ ~ 

and r\+/2 ~ Mio are fine-tuned. To be consistent 
with these scalings the EFT couplings in Eq. (l25l l must be- 
have as 



^1+ 



lo 

2n 



1 



3;r fi^Mio ' 



and 



4 



MioMhi 
(29) 



0+ 



2fi 



and 



00+ 



1,2 



(30) 



in 5 1/2, where /i now stands for the pa reduced mass. In 
the latter case the situation is nearly identical to the proton- 
proton case y_9J; up to NLO the 51/2 amplitude reads 



27r 



yu -l/flo+ -2kcH(ri) 



nuk^/2 



-l/ao+ -2kcH(j]) 

(31) 



For P-wave interactions. Coulomb can be introduced 
in a straightforward way ||3T|| . The evaluation of Feynman 
diagrams and the necessary resummation results in 



(1)2 2,0-1,2 



1 + 



1 + 



where 



2;r 

'Pi^kV4 - 

-l/ai^+ri^ky2-2kcHm(7l) 



H^'\i^)^k\\+if)H{r^), 



(32) 

(33) 
(34) 



and the variable kc - ZaZpfiagm is adapted to the pa case. 
The f 3/2 projector is given by 



Pi+{9) = 2 cos + o- ■ n sin , 



(35) 



where 9 is the angle between k and k' — the initial and final 
momenta in the center-of-mass frame, respectively — and 
n = kxk'/lkxk'l is the unit vector normal to the scattering 
plane. In this form, the T\+ amplitude requires the same 
kinematical fine-tuning discussed in Sect. 13.21 For prag- 
matical reasons, we therefore adopt the expansion around 
the resonance pole in this channel. 



4.1 P3/2 resonance pole expansion 

For practical applications, the more efficient way of ex- 
pressing the EFT amplitude in the presence of low-energy 
resonances is to perform the resonance pole expansion. 
One nice feature is that, at every order in the power-counting 
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the amplitude maintains its resonant behavior at the ex- 
pected position with the correct magnitude. It is straight- 
forward to extend the ideas outlined in Sect. 13.21 to the 
present case, with a resonance pole in the fourth quadrant 
of the complex momentum plane 1 3 1 1 . With the pole given 
by Eq. ( |28] ) one rewrites Eq. ( l32l l as 



In 



'f(k^-kl)-2kc 



H(i)(k)-Ha)(k) 



X < 



1 + 



'^(k^-kp-lkc 



H(i)(fc)-7/a)(fc) 



where 

LO NLO 

and Li, L,- are defined at the pole position from 
2kcH^^\kclkp) ^klLr + likjki Li . 



(36) 



(37) 



(38) 



The factors in front of Lr, Li ~ 0(1) reflect the behavior of 
the Coulomb function H^^\t]) around the resonance. No- 
tice that, as in the aa amplitude, Ti+ is parametrized only 
in terms of k,- and A:, at LO, with the extra parameter !Pi+ 
appearing at NLO. 

The P3/2 scattering length and effective range are given 

by 



2 



1 

fli+ 




Mf and ri+/2~Mto 



and nicely show the scalings l/ai+ - ir-ij^ 
that were expected from the numbers of Table [3] 



4.2 results 

In order to test our power-counting assumptions, we fit the 
expressions for the S 1/2 and P3/2 amplitudes to the corre- 
sponding phase shifts and compare the obtained ERE pa- 
rameters with the ones in the literature. We used the num- 
bers from Table 5 of Ref. |33i in our fit. The results are 
shown in Table H] In our fits we have assumed that all 
the points have the same statistical weight, since error bars 



Table 4. 5 1/2 and Pj/t pa ERE parameters extracted in our fits at 
LO and NLO. 



order 


ao+ (fm) 


ro+ (fm) 


LO 


7 4+8.0 




NLO 


^•°'^-0.21 


1 7+'-3 
'-0.8 



order 


ai+ (fm^) 


n+ (fm-') 


Pi+ (fm) 


LO 


-58.0+2g;55 






NLO 


-44 5+'-'' 


"■^"-0.10 


T 9 + 1.0 



Table 5. pa resonance parameters extracted from our fits at 
LO and NLO. 



P3/2 


k,- (MeV) 


ki (MeV) 


Er (MeV) 


r«/2 (MeV) 


LO 




10 3+'* 

iO.J_Q 8 


1 64+" "' 


70+° '^ 

'"-0.09 


NLO 


50.7!°;^ 


9 40+" '" 


1 66+" ""• 

^ -""-0.04 


63+°"' 

"■"-^-002 



were not provided. However, we incorporated the EFT sys- 
tematic error estimates by imposing that the amplitude can- 
not be reproduced to a precision better than {k/2m„)"^^ , 
where k is the CM momentum of the pa system, m„ our 
high-energy scale (the pion mass), and n the order in the 
EFT expansion. The maximum and minimum values of 
each ERE parameter were determined by multiplying the 
amplitude by a factor of 1 + x{k/2m„)"^\ where x is a con- 
stant that varies between -1 and 1. For the central value, 
such constant is zero. As one can see, our results are con- 
sistent with the ones shown in Table [3] at NLO. At LO the 
central values are slightly off, nevertheless consistent with 
an expansion parameter of the order of 1/3. 

In Table |5] are shown our fit results for the P3/2 res- 
onance parameters. One should stress that A:,, and A:, are 
the variables determined from the fit, from which fli+ and 
ri+ can be extracted via Eqs. ( l40t and ( l39b . The results 
are consistent with the values in Table |2l except for A:, 
which is slightly off. That is probably due to the fact that 
it is a subleading effect compared to A;,., k, ~ M^^/Mi,i, 
and to improve the agreement one needs to go one or- 
der higher. Nevertheless, the resonance energy and width 
are consistent with Ref. |32|, which used the extended R- 
matrix analysis. The latter is based on the same concept 
that the resonance is a manifestation of a pole of the am- 
plitude in the lower half of the complex energy plane. Both 
our and Ref. 132] results disagree with the traditional R- 
matrix analysis, which are nowadays considered less reli- 
able in extracting resonance properties |30|. 

In Fig. [3] we show the EFT results for 51/2 and P3/2 
phase shifts at LO and NLO. A very good convergence is 
observed, especially for the P3/2 channel. The curves using 
the ERE expressions are omitted, since they essentially fall 
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on top of the NLO curve and cannot be distinguished in a 
plot. 




Ep [MeV] 

Fig. 3. EFT results for 5 1/2 and Pj,i2 scattering phase shifts at 
LO (dotted) and NLO (sohd), compared against the partial wave 
analysis results from Amdt et al. (diamonds). 



Since we obtained a very good fit of the 5 1/2 and f 3/2 
EFT amplitudes to the corresponding phase shifts, it is use- 
ful to compare directly to observables in order to test our 
assumptions about the neglected higher-order terms. This 
is shown in Fig. |4] for the elastic differential cross-section 
measured at the laboratory scattering angle 6 - 140°. The 
sharp enhancement at low energies reflects the characteris- 
tic dominance of the Mott cross-section. One clearly sees 
that already at LO (dotted line) one has a good description 
of data, especially around the resonance peak. There is a 
slight improvement at NLO (thick solid line). The curves 
also agree well with the ERE curve (thin solid line), which 
uses the ERE formula for the 5 1/2, P^ii, and P1/2 partial 
waves and the numerical parameters from Ref. |33 1. How- 
ever, small deviations start to show up at Ep » 3.5 MeV, 
indicating an increasing importance of higher order terms 
as one goes higher in energy. The deviations, which seem 
to come mainly from the P1/2 channel, can be eliminated 
by going to higher orders in the power counting. 



5 summary 

I reviewed the programme of halo/cluster EFT suitable to 
study several interesting nuclear processes that play a role 
in nuclear astrophysics. I concentrated on two important 
aspects commonly present in those systems, namely, nar- 
row resonances and Coulomb interactions. The aa and pa 
interactions were given as examples of how setting up the 
formalism, then used as practical applications. 

The study of the aa interaction revealed a very interest- 
ing picture, where the scales of the strong interaction con- 
spire to produce a system with an almost non-relativistic 
conformal symmetry. One also observed extra amount of 
fine-tuning once electromagnetic interactions between the 
particles are turned on, up to the point of generating an 5 - 
wave scattering length that is three orders of magnitude 




2 3 4 
Ep[MeV] 

Fig. 4. EFT results at LO (dotted) and NLO (thick solid) for pa 
elastic cross-section at 6 = 140° laboratory scattering angle, com- 
pared against the partial wave analysis results from Arndt et al. 
(thin solid) and measured data point from Ref. L34J (diamonds). 



larger than the typical range of the interaction. Despite 
the struggle in understanding the subtle cancellations be- 
tween the strong and electromagnetic forces, the formalism 
is quite successful phenomenologically. In fact, thanks to 
those cancellations we were able to pin down the 5 -wave 
effective range parameters with an improved accuracy rel- 
ative to previous studies. 

We extended the formalism to include P-waves with 
resonance and Coulomb interactions when dealing with 
the pa system. We performed the expansion of the f 3/2 
amplitude around the resonance pole, which allowed us 
to extract the resonance properties directly from a fit to 
the phase shift. We derived expressions for the scattering 
length and effective range in terms of the resonance param- 
eters, that explains the scalings of the former with the low 
and high momentum scales M/o and M/„ of the theory. Our 
results at LO and NLO exhibit good convergence at a ratio 
of about 1/3, and the resonance energy and width are con- 
sistent with the ones using the extended /^-matrix analysis. 
Comparison with the differential cross-section at 140° lab- 
oratory angle reassures the consistency of the power count- 
ing, with Pi/2 contribution starting to show up only for pro- 
ton energies beyond 3.5 MeV. 

The aa and nucleon-or interactions are the basic ones 
before considering more complicated clusters of a and nu- 
cleons. An interesting example is the Hoyle state in '^C, 
which plays a key role in the triple-alpha reaction respon- 
sible for the formation of heavy elements. A model-study 
with this state in mind was developed in Ref.|35|, where 
a perturbative treatment of the Coulomb interaction was 
proposed. This idea might be useful to handle the techni- 
cal difficulties involving three charged particles. 
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